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A STUDY OF THE HARTREE-FOCK MODEL SPACE FOR 


LIGHT DEFORMED NUCLEI 

by William F. Ford, Richard C. Braley, 
and J. Bar-Touv* 

Lewis Research Center 
SUMMARY 

The deformed Hartree-Fock method is used to study the effect on various nuclear 
properties of altering the truncated basis space used in structure calculations. After a 
brief discussion of the Hartree-Fock model and the means by which nuclear properties 
are obtained from the projected physical states, results are presented for the nuclei 
neon-20, magnesium-24, silicon-28, and sulfur-32. Intrinsic spectra, nuclear radii, 
electromagnetic transition rates, electron scattering form factors, and inelastic proton 
scattering cross sections are investigated and wherever possible compared with experi- 
ment. Calculations are performed with both Wood-Saxon and harmonic oscillator basis 
functions; 6 Wood-Saxon functions are used, and up to 15 harmonic oscillator functions 
with various values of the oscillator length. Results show that changes in the basis 
space produce rather dramatic changes in the nuclear properties, indicating that 
Hartree-Fock calculations of the type presented here are being performed in an inade- 
quate basis space. 


INTRODUCTION 

The shell model and many of its ensuing developments have established that a mi- 
croscopic description of the nucleus, in terms of the motions of individual nucleons, 
can be successfully used to predict energies, spins, and parities of the low-lying nuclear 
states (refs. 1 to 3). While such descriptions are generally designed to study level sys- 
tematics, the wave functions so obtained have occasionally been subjected to more rigor- 
ous tests, such as prediction of electromagnetic transition rates (refs. 4 to 6), electron 
scattering form factors (refs. 7 and 8), and inelastic proton scattering cross sections 
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(refs. 9 to 11). Most of these properties are quite sensitive to the behavior of the 
wave function in the vicinity of the nuclear surface, whereas energy level systematics 
are thought to be more dependent on details within the nuclear interior. Calculations 
of surface-dependent properties have been, for the most part, less successful in re- 
producing experimental data. 

In recent years the number of such microscopic investigations has steadily in- 
creased. These studies differ from one another mainly in the type of two-body interac- 
tion employed and in the size of the model space. The interaction may be designed to 
fit the nucleon-nucleon scattering data (ref. 12), may be of the renormalized or "ef- 
fective” type (ref. 13), or may be obtained from an approximate solution of the Bethe- 
Goldstone equation (refs. 14 and 15). The model space nearly always employs harmonic 
oscillator wave functions, but there is little uniformity in the choice of oscillator length. 
Enough functions are included to allow for configuration mixing or deformation of orbits; 
for the heavier nuclei, calculations are often kept within reasonable bounds by presum- 
ing the existence of an inert, closed-shell core. 

Efforts to extend the early shell model calculations have been devoted primarily to 
spherical vibrational nuclei, with correlations between nucleons moving in an average 
field being introduced via pairing (ref. 2). Excited states of these nuclei can be ob- 
tained by using some form of the random-phase approximation (refs. 2 and 3). 

At the same time there has been significant progress in our understanding of de- 
formed rotational nuclei. Here the chief means of investigation has been the Hartree- 
Fock (HF) method, which has received its greatest attention in the study of light de- 
formed nuclei. One of the major advantages of working with light nuclei is that the inert 
closed-shell core assumption can be dropped and the problem is still tractable with re- 
spect to projection of states of good angular momentum from the intrinsic HF wave 
functions. 

Recent calculations by several groups indicate that the HF method (or HF with 
slight variations) provides a good description of low-lying states of the light deformed 
nuclei (refs. 1, 16 to 18). In most of the nuclei studied it was found that the correct 
level ordering was reproduced, with the level spacings in good agreement with experi- 
ment (refs. 5, 16, 17, and 19). On the other hand, E2 rates and inelastic proton cross 
sections obtained from projected HF states indicate that these functions do not have the 
proper behavior needed to predict correctly the magnitudes of quantities associated with 
inelastic transitions. If the underlying assumption of the HF method is correct, namely, 
that excited levels within a band are generated from the same intrinsic state, then the 
source of the difficulty may well be the asymptotic form of the intrinsic HF state itself. 

Most HF calculations have made use of basis functions with harmonic oscillator 
radial dependence, in spite of their rapid fall- off at large distances. Of course, if a 
sufficiently large number of such states were included, the correct asymptotic behavior 
could be produced. Indeed, recent studies by several groups indicate that calculated 
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nuclear properties improve - relative to experiment - as the size of the model space 
increases. Bassichis, Pohl, and Kerman (ref. 20) have studied the intrinsic HF state 
using both four and six major shells. They found that the quadrupole moment is quite 
sensitive to extension of the space, while the binding energy, rms radius, and kinetic 
energy per particle apparently are not. Gunye (ref. 5) examined projected energies for 
several nuclei in the lower half of the 2s-ld shell and found the spectra improved con- 
siderably as the space was enlarged to include four major shells. Preliminary inves- 
tigations by the present authors (refs. 21 and 22) revealed that E2 rates and inelastic 
proton cross sections obtained from projected HF wave functions are substantially af- 
fected by the number and kind of basis functions used. These findings provide ample 
justification for further study of the effect on HF calculations of variations in the model 
space, particularly with regard to inelastic transitions induced by scattering or interac- 
tion with the electromagnetic field. 

Microscopic descriptions of inelastic nucleon scattering have received considerable 
attention in recent years (refs. 9 to 11) largely as a result of advances in the theory of 
nuclear structure. Most applications have been to nuclei in the collective vibrational or 
rotational regions where a large number of structure calculations have been made. The 
distorted wave Born approximation (DWBA) is generally used to describe the scattering 
process, and the structure problem is solved with an inert core. Exceptions include 
coupled channel calculations by Glendenning (ref. 11) (but with an inert core) and DWBA 
studies by the present authors (ref. 9) in which core effects are explicitly treated. 

There are indications that DWBA may be adequate for low-lying states in most of the 
light deformed nuclei, but that coupled channels will be required for states with higher 
excitation energy. 

The scattering of high-energy electrons from nuclei also provides a powerful tool 
for nuclear structure studies. As a means by which nuclear models may be tested, 
electron scattering has the attractive feature that the electron-nucleus interaction is 
well understood. This is in contrast to the situation in nucleon scattering, where uncer- 
tainties in the nuclear model are further complicated by lack of knowledge of the 
nucleon-nucleus force. 

The analysis of high-energy electron scattering is carried out most easily in the 
first Born approximation; however, it is well known that this method is inadequate for 
large momentum transfers and that zeros occur at the diffraction minima. Baker 
(ref. 23) has proposed an approximation for high-energy electron scattering based on 
the eikonal approach developed by Glauber (ref. 24). Baker's high-energy approxima- 
tion fills in the zeros of the diffraction minima and agrees well with the results of phase 
shift analyses. To our knowledge, however, this technique has only been tested for the 
simplest nuclear models. 

There is a large amount of experimental data available for the 2s- Id shell nuclei. 

As a result of this, together with evidence that many of these nuclei have deformed 
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equilibrium shapes, the 2s-ld shell has become a testing ground for the deformed HF 

20 24 

method. Of the nuclei in this region, neon-20 ( Ne), magnesium-24 ( Mg), silicon-28 

OO Q O 

( Si), and sulfur-32 ( S) have received the greatest attention because of their conven- 
ient symmetries (even-even, N = Z). While it is generally conceded that the simple HF 

model (i. e. , a single Slater determinant for the intrinsic state) provides a fairly good 
20 24 2ft 

description of Ne and Mg, the simple model appears to be inadequate for Si and 

OO 

S. Accordingly there have been several attempts to explain their structure in terms of 
slight modifications of the HF method (refs. 16, 17, and 25). Among these, Bar-Touv 
and Goswami have proposed a model of "inverted coexistence" which seems to be prom- 
ising (ref. 25). It remains now to make microscopic studies of these nuclei with ex- 
tended model spaces. 

The purpose of this work is to examine, as consistently as possible, the effect of 
alterations in the model space on nuclear properties predicted by the Hartree-Fock 
method. Wave functions are obtained by projecting states of good angular momentum 
from intrinsic HF states. Several nuclear properties are examined for model spaces of 
different size and with Wood- Saxon radial dependence as well as the usual harmonic 
oscillator behavior. 

For the comparison of Wood-Saxon and harmonic oscillator cases only six basis 
functions are used, since all the other Wood-Saxon eigenfunctions are unbound. In this 
case the severe truncation makes it necessary to use an "effective" nucleon-nucleon 
interaction, and the Volkov force (ref. 26) was chosen. For the comparison of model- 
space size, sets of 6, 10, and 15 harmonic oscillator functions are used with varying 
oscillator lengths. Here one of the "realistic" nucleon-nucleon interactions, namely, 
that devised some years ago by Tabakin (ref. 27), is used. 

In the next section, a brief review of the HF method is given, and a summary of the 
technique used to calculate nuclear properties from the projected wave functions. De- 
tails of the calculation are outlined, the results are discussed, and a summary and some 
conclusions are presented in later sections. 


THEORY 
Nuclear States 

The low-energy spectra of some nuclei remind one of the spectrum of a rigid rota- 
tor, that is, the energy levels are proportional to J(J + 1). This observation led to the 
early descriptions of deformed nuclei in terms of a rotating intrinsic state (ref. 28). 
Thus, a natural starting point for the deformed HF method was provided. In the HF pic- 
ture, the intrinsic state is described in terms of nucleons moving in a deformed field 
with respect to body-fixed axes. Physical states of the system may be thought of in 
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terms of various rotations of the intrinsic state. Formally, the physical states may be 
obtained from the intrinsic state by the well-known method of angular momentum pro- 
jection (ref. 29). 

Intrinsic states. - In the HF method the intrinsic nuclear wave function is approxi- 
mated by a single Slater determinant 

$ K = (A! T 1 / 2 det{<p x (x n )} (1) 

whose orbitals <p^(x n ) are determined by minimizing the expectation value of the intrin- 
sic nuclear Hamiltonian. Solution of this variational problem leads to the HF equations 

h|x>=ejx> (2a) 


where 


<i|h|j> = < i 1 1 1 j) + E < iX | V | jX) (2b) 

A 

The i and j refer to single-particle (s.p. ) states, and the sum on X runs overall 
occupied orbits. In equation (2b), t represents the s. p. kinetic energy operator, and 
the tilde indicates that the two-body matrix elements are antisymmetrized. 

To solve the HF equations, each of the orbits is written in the form of a (truncated) 
expansion 


|*> = £c*|j> (3) 

j J 

in which the summation over j is partly restricted by the symmetries of the nuclear 
system. The are determined self-consistently from equations (2), and the HF 
energy is then found from the expression 

e hf 4 2 [ < " |t|t ‘ > +e J (4) 

Ph ysical states. - The intrinsic state obtained from the solution of the HF 
equations is not an eigenfunction of total angular momentum. However, there is a well- 
known prescription for obtaining such eigenfunctions, namely (refs. 1 and 29), 
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( 5 ) 


^ JM ( 1, • • • > A) = ^L±i /dS2 D^ k (S2)*R(S2)$ k (1, . . . , A) 

8jt 2 

where is the rotation matrix, R is the rotation operator, and M and K are 

the projections of the total angular momentum J onto the space- and body-fixed axes, 
respectively. When the restriction to axial symmetry is made on the HF state, equa- 
tion (5) simplifies to 


/ ~idJ 

sin 9 dd d^ K (0)e y # K (1, A) (6) 

^ U 

where is the reduced rotation matrix. In all of the analysis which follows the 

nuclei to be considered have even N and even Z, and only states within the ground state 
rotational band are considered (i. e. , K = 0). Consideration of the symmetries of <t>^ 
enables one to show that, for even-even nuclei, the J values of states in the ground 
state band are all even. 


Nuclear Properties 

The prediction of nuclear properties requires evaluation of reduced matrix elements 
of certain operators between initial and final states of the nucleus. Most of these oper- 
ators can be expressed as a sum of single-particle operators (an exception, of course, 
being the nuclear Hamiltonian). When this is the case, it is possible to evaluate the re- 
duced matrix elements in such a way that the integration over the coordinates of the 
single-particle operator is the last to be performed. The previous integrations over the 
other (A-l) coordinates thus yield a quantity which depends only on the initial and final 
states of the nucleus and the nature of the transition involved. This leads directly to the 
concept of an inelastic transition density p^gj(r), which essentially gives, for the trans- 
ition i — f, the probability for absorption of angular momentum J = L + S by the last 
nucleon as a function of its distance from the nuclear origin. This very useful concept 
has been discussed in detail in a previous publication (ref. 30). 

In the cases considered here, each single-particle operator O may be expanded in 
terms of spin-angle tensors with coefficients co L gj(r) which give the radial dependence 
of the operator. Reduced matrix elements of the many-body operator may then be ex- 
pressed as 
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^LSJ (r)co LSJ (r)r2 dr 


(7) 



Z) fi(n) 

n=l 



The expansion of PLSj( r ) * n terms of a biorthogonal set of (Gauss- Laguerre) functions 
is particularly convenient from a calculational point of view since it enables one to re- 
express equation (7) as 


<J f IR|Ji> = Zt p LSJn f ^ dr (8) 

LS n=0 * / o 

For most operators the integral over r can be carried out analytically, and the sum 
over n has been found to converge rapidly. This is accomplished by suitable choice of 
the adjustable parameters a and /3 appearing in the Gauss- Laguerre functions. 

Static nuclear properties. - The initial and final states of the nucleus are the same 
in this case. For ground state properties of even-even nuclei, the calculation of the ma- 
trix elements is particularly simple since = 0 and L = S = J = 0. One of the 

ground state properties which we require to be in agreement with experiment is the rms 
radius. In this case, the operator is 


R 


2 



n 


(9) 


and equation (8) becomes 


<i||R 2 ||i> 





p 000n 


where, for convenience, we have introduced 


~if = ( a \ n ^if 
p LSJn ~ P LSJn 


( 10 ) 


(U) 


The electric quadrupole moment provides information about the shape of the nucleus. 
To obtain it we evaluate the reduced matrix element of 
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( 12 ) 


Qo = e 


(16 * 


. 1/2 


Z) r ? Y 2( f i> 

i=l 


1 + T t 


th 

where r n is the z-component of the isotopic spin operator for the n n nucleon. The 
result is 

<Jfl|Q 2 ll J f> - e(f) V V 2 (l - »r 7/2 L~4% n (13) 

\ o ) n=0 

Only the proton part of the transition density is involved here, as is indicated by the p 
superscript on p. 

As pointed out in the INTRODUCTION, the elastic scattering of high-energy elec- 
trons is described more accurately by the Baker approximation than by the first Born 
approximation. The form factor for electron scattering is given by 

da 

I F(q) | 2 = -y— r — (14) 

(s)» 

in which (do/dfi)j^ is the Mott cross section and 

| M | 2 


where 



(15) 


as shown by Baker. When microscopic wave functions are used to represent the target, 
the evaluation of equation (15) is rather involved; it is presented in the appendix. 

Dynamic nuclear properties. - The prediction of quantities which involve transitions 
to excited states provides a more sensitive test of the validity of any nuclear model. Of 
particular interest to us shall be electromagnetic transition rates and inelastic particle 
scattering. For the EX transition rate from i to f we obtain (ref. 30) 


8 


2 


(16) 


B(EX; i - f) = e 


2 2Jf + 1 

2J. + 1 


r * (1 . CT) -M3/ 2 )£~«P n 


For the inelastic scattering of high-energy electrons we evaluate the Born approximation 
to the electron scattering form factor, obtaining an expression of the form 


F(q) 


2 


2 J f + lVl V* 

2 J i + 1 ^| L n=0 


^LoLn^nL^ 1 " a '> 


2 


r'q) 


(17) 


Details of the derivation may be found in reference 30. 

For the inelastic scattering of protons it is necessary, as mentioned in the INTRO- 
DUCTION, to make some assumptions concerning the nature of the nucleon-nucleus 
force responsible for the transition. The standard assumption, and the one adopted 
here, is that this interaction may be represented as a sum of two-body M effective" po- 
tentials: 


A) = £ v(x x n ) (18) 

n=l 

In order to calculate the scattering amplitude, the reduced matrix elements of 

V(x , A) must be obtained as functions of the projectile coordinates x , and then ex- 
o if ° 

panded in spin-angle tensors; the radial coefficients f;Lgj( r 0 ) are known as th e nucleon 

scattering form factors (ref. 9). Because of the above assumption, V is a one-body 

operator in the coordinates of the target nucleons and the methods of equations (7) and (8) 

apply. The resulting form factors are expressible as 

f LSj( r o) = f LSJN*^NL^ 0, ?r o) 

For details the reader is again referred to reference 30. 


DETAILS OF THE CALCULATIONS 

The HF variational problem was solved subject to the condition that the intrinsic 
state possess axial symmetry and fourfold degeneracy (ref. 1). Each orbit is therefore 
specified by its z-components of spin and isospin and further specified by its projections 
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onto basis functions labeled by the quantum numbers (nZ j). The model spaces to be used 
in this work are listed in table I. 

The final choice of oscillator length b = for the harmonic oscillator func- 

tions is based on the requirement that the rms radius be well reproduced for each nuc- 
leus, that is, the rms radius serves as a constraint in the variational problem. How- 
ever, the effect (on other quantities) of variation of the oscillator length is also investi- 
gated. 

The parameters of the potential which produced the Wood-Saxon functions were 
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chosen so as to yield approximately the 0 experimental single-particle energies 
(ref. 4). This potential has the form 


V(r) = V^r) - aV Q 



2 


£ ± P(r) 

r dr 


with 


p(r) = 



-1 


and V Q = -56. 5 MeV, R = 3. 15 femtometers, a = 0. 65 femtometer, and a = 20. There 
are states generated by this potential other than those shown in table I, but these lie in 
the continuum. 

The two-body matrix elements have been calculated using the Tabakin interaction 
for most of the cases in which harmonic oscillator functions are used. This is a smooth, 
nonlocal separable potential which has been used successfully in several shell model and 
HF studies (refs. 12 and 20). The general form of the potential is 


rM 


' h » l (r)l v ’ <r ' )} Y “ CrY C 

ll' 


details may be found in reference 12. 

For calculations involving the comparison of Wood-Saxon and harmonic oscillator 
functions, one of the "effective” nucleon-nucleon interactions was used, as mentioned 
in the INTRODUCTION. We have chosen a Volkov force (ref. 26) which yields fairly 
good saturation properties for the p- shell nuclei: 



where 


f(|r t - rj) = - 78 e 


rr 2 

r l r 2 


)/2. 25 


+ 82. 5 e 


-( 


- - 2 
r l -r 2 


)/0. 64 


with 0 = 0. 1875, /3 = -0. 0775, y = 0. 2025, and 6 = -0. 1775. Another reason for the 
choice of this potential is that it is local, which is particularly convenient for calculat- 
ing two-body matrix elements with Wood-Saxon wave functions using the technique of 
reference 4. 

Transition amplitudes for inelastic proton scattering are obtained in the DWBA, 
neglecting exchange effects between projectile and target nucleons. The following opti- 
cal potential is used to generate the distorted waves: 

U(r) = V c (r) - V Q Py(r) - ^ W s p 8 (r) - V LS (^-) 2 ^_p LS (r) 

a _ \ rn c / 3,-r c?r 


where V is the Coulomb term, V the strength of the real potential, W 0 the strength 
CO s 

of the surface absorption, and V LS the spin- orbit strength. The quantity Py(r) is the 
usual Wood-Saxon form factor, and pg and p^g are derivatives of functions similar to 
Py but with appropriate radius and diffuseness parameters. 

Two different forms of the interaction between projectile and target nucleon are 
considered, the first due to Glendenning and Veneroni (ref. 10): 


-(r 12 /1.85) 2 

V( r i 2 ) - -52 e (P T E + 0-6 PgjO 


and the second (with a considerably shorter range) to Thompson and Tang (ref. 31): 


-(r 12 /1.47) 2 

V(r 12 ) = -72.98 e (P TE + 0. 632 Pg E ) 


RESULTS 

Harmonic Oscillator Basis 

In the usual structure calculation which employs harmonic oscillator basis functions, 
the number of functions to be included and the value of the oscillator length are ’’free 
parameters. ” Although the effect of enlarging the model space has been investigated in 
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some detail, the effect of variations in the oscillator length has not been studied exten- 
sively. When a truncated basis is employed, however, the results will depend on the 
oscillator length used, and the nature of this dependence should be examined. 

In this section we present results indicating the effect on various nuclear properties 
of varying both the number of basis functions and the value of the oscillator length. In 
the final analysis, the oscillator length is adjusted so that the rms radius for each nu- 
cleus, calculated in the largest model space, is in agreement with experiment. 

Neon-20 . - We begin by presenting results for fixed oscillator length. In table n(a) 
are compared the intrinsic HF states and single-particle energies for the three model 
spaces defined in table I, using an oscillator length of 1. 88 femtometers. In table 11(b) 
the HF energies and gaps (between last occupied and first unoccupied orbitals) are also 
compared. The level ordering is preserved in all three cases, and the HF energy de- 
creases (binding energy increases) as the model space is enlarged, as has been observed 
previously (ref. 20). The energy gap A increases as well, indicating that the HF so- 
lution becomes more stable as the model space is enlarged. 

We now turn to the physical nuclear properties which are obtained by projecting 

states of good angular momentum from the intrinsic HF state. The projected energy 
20 

spectrum for Ne has been shown to improve as the model spacing is enlarged (ref. 5). 
However, other important nuclear properties have gone untested in this regard. 

In table III we show the effect on the projected rms radius and E2 transition rate of 
variations in both the size of the model space and the value of the oscillator length. Ta- 
ble in(a) gives the results for fixed oscillator length. As the space is enlarged, the E2 
rate increases, although even in the largest space it is far below the experimental value. 
On the other hand, the nucleus contracts as the number of basis functions is increased 
and eventually has a radius smaller than the experimental value. (Experimentally, 

2 . , o 4 2 

<R ) = 2. 79 fm, and B(E2; 0 + - 2 + ) = 286 e^ • fm“ The error on (IT) is 

probably of the order of 5 percent. Experimental radii were obtained from R. de 

Swiniarski, and E2 rates appear in ref. 32. ) 

Table in(b) gives the results for a fixed number of basis functions, showing that as 
one might expect both the rms radius and the E2 rate increase with oscillator length. 

This is due to the displacement toward larger radius of the maxima of the radial func- 
tions as b is increased. 

Finally, table HI(c) shows the results when the rms radius is essentially fixed, 
serving as a constraint on the variational calculation. It should be noted that in none of 
these cases does the value of b minimize the HF energy, considered as a function of b. 

A typical example is shown in figure 1, where the HF energy is plotted as a function of 
oscillator length (for space 2). A smooth curve is drawn through the points indicating a 
minimum HF energy at b = 1. 83 femtometers, whereas table in indicates that the cor- 
rect rms radius is achieved at b = 1. 88 femtometers. 
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As an additional test of the wave function, the E2 rates between other states within 
the ground state band are predicted and compared with experiment in table IV. The E2 
rate for 4 + — 6 + is underestimated in much the same fashion as the 0 + — 2 + rate; how- 
ever, the 2 + — 4 + rate is slightly beyond the experimental limit. It is surprising indeed 
to find a predicted E2 rate which overestimates the experimental result; in fact, in 
view of the experimental uncertainty in the rms radius, it is not unreasonable to claim 
agreement with experiment for this transition. Of course, the agreement may be for- 
tuitous; but if not, the result suggests that perhaps each of the physical states should be 
generated from a different intrinsic state. Recent investigations by Castel and Parikh 
(ref. 16) indicate that this may be a better procedure than the standard approach, in 
which the idea of a band generated from a single Slater determinant is invoked. 

Unfortunately, there are no experimental data available, at the present time, from 

20 

which elastic electron scattering form factors for Ne may be extracted. However, 
since the method developed by Baker has not been used previously to investigate micro- 
scopic nuclear wave functions, it seems worthwhile to display some of the changes 
which result from variations in the nuclear model and in the parameters which describe 
the scattering. In figure 2 the form factors are compared (for space 3 with b= 1. 93 fm) 
in the Baker approximation and in the Born approximation. The more sophisticated 
method fills in the diffraction minimum and shifts the second maximum to lower q. 

Also shown is the result which occurs if correction is not made for the finite size of the 
proton. The effects of enlarging the model space are shown in figure 3; the trend is to 
shift the first minimum to larger q and to increase the magnitude for smaller q 
values. 

We turn now to the case of inelastic proton scattering, recalling that the structure 
of the target nucleus as well as the radial dependence and mixture of the two-body force 
manifest themselves in the shape and size of the nuclear form factor. In figure 4 the 
scalar form factor fj^M is presented for the 0 + — 2 + transition, using the 
Glendenning-Veneroni interaction. Although the shape of the form factor does not change 
noticeably, its magnitude increases with the number of basis functions. The vector form 
factor ^ as rou S h ly the same shape as the scalar but is usually smaller by about 

two orders of magnitude. 

The cross sections obtained using these form factors are shown in figure 5. As the 
model space is enlarged, the cross section nearly doubles, although its shape is not 
changed significantly. The angular distribution is given reasonably well; however, it 
underestimates the experimental data (of de Swiniarski) by at least a factor of 2 in the 
forward direction. 

Since there are uncertainties in the interaction between the projectile and the bound 
nucleons, these calculations are repeated for the Tang force. As can be seen from the 
comparison in figure 6, the major difference in the form factors occurs for radii larger 
than 4 femtometers. The effect which this has on the cross section is shown in figure 7: 
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the magnitude of dcr/dfi using the Glendenning-Veneroni force is generally larger than 
that using the Tang force, despite the fact that the latter is about 40 percent stronger. 
This occurs because of the short range of the Tang force, which makes the form factor 
smaller in the important surface region. 

24 

Magnesium-24 . - The early HF studies of Mg by Bar-Touv and Kelson indicated 
that the nucleus is asymmetric, because of the larger binding energy and gap possessed 
by the asymmetric solution. However, comparison of projected spectra for the two 
shapes (ref. 19), together with recent studies in which pairing effects are included 
(ref. 33), indicates that its intrinsic structure is still in doubt. Furthermore, Bar- 
Touv has observed that the symmetric and asymmetric states have a large overlap. For 
purposes of this study, therefore, there is little reason to prefer one solution to the 
other; the symmetric solution is actually chosen for convenience in projecting states of 
good angular momentum. 

The intrinsic HF states and single-particle energies are compared in table V(a). 

24 

(For brevity, results in the smallest space will not be presented for Mg on the re- 
maining nuclei.) The HF energies and gaps are shown in table V(b), and the variation 

of the HF energy with oscillator length in figure 8. The general trends are clearly the 
20 

same as for Ne, and we omit much of the detail; results for the ultimate choice of 
oscillator length are given in table V(c). 

Some of the electromagnetic transition rates predicted are presented in table VI; 

20 

the 0 + — 2 + rate is underestimated by about 25 percent, but as in the case of Ne, the 
2 + — 4 + is overestimated. It appears that the wave functions for both of these nuclei 
have similar inadequacies. 

The electron scattering form factor is compared with unpublished experimental data 
of H. Hultzsch in figure 9. The predicted result is quite good for q-values up to 1. 2 
femtometer - ! but does not correctly predict the first minimum nor does it rise suffi- 
ciently at the secondary maximum. It is not clear whether this is due to the incorrect 
shape of the Mg wave function in the nuclear interior or some inadequacy in the method 
for computing the form factor. 

The DWBA prediction for inelastic proton scattering to the first 2 + level is compared 
with experiment (ref. 34) in figure 10. The general quality of the predicted cross section 
is the same as for Ne: the angular distribution is fairly good but drops off a little too 
rapidly at larger angles, and the prediction underestimates the data by roughly a factor 
of 2. 

Silicon-28. - As mentioned in the INTRODUCTION, there are several reasons for 
suspecting that Si is too complicated for description by the simple HF picture, and 
that one must resort to the more involved methods such as "inverted coexistence" or the 
inclusion of various other types of correlations. The projected energy spectrum yielded 
by the HF model is far too compressed, for one thing, and the 0 + — 2 + transition rate is 
badly underestimated. 
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On the other hand, all spectra and E2 rates calculated to date have made use of 

limited model spaces. Enlarging the model space is known to have a significant effect 

on these quantities, so there is still reason to study this nucleus in the simple HF pic- 

28 

ture. In the discussion of Si, we will first examine the results obtained from the HF 

model, and then present some results based on the assumption of inverted coexistence. 

The oblate solution of Si is the lowest of the three (the other two being prolate, 

with Ejjj, = -58. 30 MeV, and spherical, with Ejjp, = -45. 29 MeV, for the largest space 

with b = 2. 09 fm), and unless otherwise specified, will be the one for which results are 

given. In tables VII(a) and (b) the HF intrinsic states, HF energies, and energy gaps 

20 

are compared for spaces 2 and 3. The general trends are similar to those for Ne and 
24 

Mg. Results obtained when the oscillator length is varied are also quite similar and 
will be omitted for ^Si and ^S. 

OQ 

In table VII(c) one finds striking evidence, however, of the difference between Si 
and the preceding nuclei. The ground state E2 rate is overestimated for both spaces, 
and decreases as the model space is enlarged. It is possible, of course, that this rep- 
resents convergence toward the experimental value from above, much as the results for 
20 24 

Ne and Mg seem to indicate convergence from below; however, the discrepancy be- 
tween the measured 2 + — 4 + transition rate (60. 1 e^fm^, e = 1. 6x10”^ C) and the pre- 
dieted value (193 e fm ) is far greater than found for the previous nuclei. The electron 
scattering form factor is shown in figure 11, and the proton scattering cross section in 
figure 12. Although the quality of the former is about the same as found for Mg, the 

drop in the proton cross section for angles smaller than 40° is discouraging. The over- 

28 

all impression is that Si is not well described by the simple HF picture. 

As mentioned earlier, there have been proposals for modifying the HF treatment of 
^Si. None of these has been tested using a detailed description in which the core orbi- 
tals are varied and an extended basis is used. The model of inverted coexistence (IC) 
proposed by Bar-Touv and Goswami (ref. 25) will be tested here. In this model, the 
0 j ground state and 0^ excited state wave functions are taken to be linear combinations 
of the spherical and deformed (oblate) HF solutions: 

\}y 1 (0 + ) = tv4> g + jS$ D 
i// 2 (0 + ) = cy’$ s + |3’$ d 

The values of a and /3 are taken from the results of Bar-Touv and Goswami, and 
and /3' are determined from the requirements of normalization and orthogonality. 

In table vm the predictions based on the IC model are compared with those pre- 
dicted by the simple HF model and with experiment. With the IC model all three quanti- 
ties are brought within the experimental limits. Unfortunately, the IC model has no 
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effect on the E2 rate for the 2 + - 4 + transition and negligible effect on the electron 

scattering form factor and inelastic proton cross section. (It should be noted that the 

proton scattering cross section might be changed considerably in a coupled channel 

calculation because of the strong coupling between the 2 + and excited 0 + states. ) 

no 

Sulfur-32. - Generally the comments made previously with regard to Si also ap- 
ply to 32 S. The only specific exception to the introductory remarks is that early HF 
studies for this nucleus indicated that it might be asymmetric (ref. 35). 

For 32 S, the prolate solution yields the lowest HF energy. Tables IX(a) and (b) 
compare the HF intrinsic states, HF energies, and energy gaps for this solution in 
spaces 2 and 3, with results very similar to those observed for Si. The radius and 
E2 rate are compared in table IX(c), again with similar findings; the result for the ob- 
late solution is also included because it is within the experimental limit (this is prob- 
ably fortuitous). The qualitative features of the electron scattering form factor (fig. 13) 
and the inelastic proton scattering cross section (fig. 14) are much like those found for 
28 Si. 

32 

An inverted coexistence calculation was not carried out for S because of the ac- 
cidental orthogonality of the prolate and spherical intrinsic states. (Our method for 
obtaining inelastic transition densities requires inversion of the matrix whose elements 
are the overlap integrals of the initial- and final-state orbits; because of the accidental 
orthogonality, the determinant of this matrix vanishes for the cross term. ) We did 

carry out such a calculation for the oblate and spherical states, however, and found it 

28 

yielded results analogous to those for a prolate-spherical Si calculation. We suspect, 

32 28 

therefore, that the IC model would probably fare no better for S than it did for Si. 


Wood-Saxon Basis 

In this section nuclear properties calculated using a Wood-Saxon basis are com- 
pared to those obtained with a corresponding harmonic oscillator basis. The choice of 
oscillator length (b = 1. 924 fm) is based on previous investigations by two of the present 
authors (ref. 4), in which it was found that this value gave the best agreement between 

nuclear properties calculated using the two different bases when the core was inert. 

20 

The intrinsic HF states and single-particle energies for Ne are compared in ta- 
ble X. The level ordering is the same for both bases, with the Wood-Saxon orbital 

energies slightly larger except for the last orbit. Similar results (not shown) are found 

32 

for these intrinsic properties of the other three nuclei, except that in S the ordering 
of the last three orbitals is different (they are nearly degenerate in the harmonic oscil- 
lator basis). 

Table XI compares some other properties for the two bases, nucleus by nucleus. 
The WS prediction for the rms radius is always smaller than the HO prediction, and, 
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except for 32 S, closer to experiment. Since the Wood- Saxon well parameters were 

17 

chosen so as to reproduce the 0 single-particle energies, it is perhaps not surprising 
that the WS prediction should underestimate the nuclear size for heavier elements. The 
E2 rates are roughly comparable for all four nuclei. The HF energies (larger here be- 
cause of the Volkov force) demonstrate an earlier (ref. 4) finding: because of the nor- 
malization to unity, the longer tails on the Wood- Saxon basis functions force them to be 
smaller within the nuclear interior, resulting in less binding energy. 

A drawback of any study in which the well parameters do not change is illustrated 

20 

in figures 15 and 16. In the former, the WS electron scattering form factor for Ne 

is seen to be very similar to that obtained previously in the largest model space, while 

32 

the HO form factor is quite poor; in figure 16, however, the reverse is true for S. 
Figure 17 demonstrates that the cross section for inelastic scattering of protons by 

20 

Ne is not significantly altered by the change in basis, except for a slight decrease in 

magnitude; the same holds true for the other three nuclei (not shown). Once again, it 

seems to be the combined action of normalization and the longer tail which is responsible 

20 

for the decrease; the effect can be seen quite clearly in the Ne form factors shown in 
figure 18. 


CONCLUDING REMARKS 

The purpose of this work has been to investigate the effects of alterations in the 
model space used for Hartree-Fock calculations. A number of different nuclear proper- 
ties have been examined; these were obtained from microscopic wave functions of good 
angular momentum projected from the intrinsic HF states. The intrinsic states were 

assumed to have axial symmetry and fourfold degeneracy, but no inert core assumption 

20 24 28 32 

was made. The 2s- Id shell nuclei Ne, Mg, Si, and S were chosen for this 
study because of their convenient symmetry properties and because of the large amount 
of experimental data available for them. 

The general results indicate that the predicted nuclear properties vary significantly 
as the basis space is enlarged. These properties are likewise sensitive to the shape and 
extent of the radial basis functions. It was observed, for example, that if the correct 
nuclear radius is to be obtained, the oscillator length must vary significantly from nu- 
cleus to nucleus; similarly, the use of a (limited) Wood-Saxon basis requires different 
well parameters for different nuclei. 

Even when a fairly large space was used in the structure study, it was found that, 

20 24 

although Ne and Mg were described reasonably well by the simple HF picture, this 

28 32 

was not true of the closed subshell nuclei Si and S. A test of the model of inverted 
coexistence for these nuclei seems to yield improved results for the ground state 0 + — 2 + 
transition rate and gives a good value for the excited state 0 + — 2 + transition rate; how- 
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ever, many other properties were unaffected by this modification, notably the other E2 
rates, the electron scattering form factors, and the inelastic proton scattering cross 
sections. 

Generally speaking, the electron scattering form factors predicted using the Baker 
approximation were unsatisfactory when compared with experiment. It is difficult to 
say whether this is due to poor nuclear wave functions or to inadequacy in the scattering 
approximation. 

Some insight into this problem may be gained by comparing the nuclear densities 
(fig. 19) with the form factors (fig. 20) for the four nuclei. One expects in general to 
find a correlation between the smaller and larger portions of these curves, and vice 
versa, simply from the nature of the Fourier transformation which they represent. The 
slope of the form factor is directly proportional to R as q — 0, for instance, and the 
frequency of the form factor minima are closely associated with lack of smoothness in 
the density for small r. 

But these curves demonstrate once again a correlation which is generally ignored, 
the correlation between asymptotic behavior and interior behavior resulting from the 
normalization condition. This illustrated rather dramatically in figure 19, where the 
reciprocal of the rms radius for each nucleus is marked along the scale at the right, to 
be compared with the nuclear density at the origin, which can be read from the scale at 
the left. The behavior of the density in the vicinity of the nuclear surface is evidently 
strongly correlated with its behavior at the origin. 

Our major conclusion, of course, is that the Hartree-Fock model gives unreliable 
results if used in a model space as small as those examined here. This is unfortunate, 
for the simplicity of the scheme is lost as the model space is enlarged, to say nothing 
of the computing time. 

On the other hand, we must point out other moderating conclusions. It is clear, 
for example, that the asymptotic part of the basis function plays an important role, and 
it may well be that the increase in the number of basis functions really represents an 
attempt to reproduce this asymptotic behavior. Therefore, no definitive statements 
about the size of the model space can be made until a study of this nature has been made 
with basis functions more appropriate to the problem. 

Another question may be raised concerning the use of a ,, realistic ,, two-nucleon 
force rather than an "effective” one. With no evidence to the contrary to guide us, we 
made the assumption that the general trends would be the same in either case. The most 
accurate procedure would be to obtain solutions of the Bethe-Goldstone equation for each 
new basis set and then proceed to the HF problem; but such a task would be prohibitively 
difficult. 

An alternative approach has been suggested recently (ref. 36). It is possible to 
develop the theory so that the operators contain the effect of the excluded space. This 
results, of course, in an effective interaction, but more importantly, it also results in 
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renormalized operators for all the observables studied here. At present the theory is 
in its preliminary stages but certainly should be considered in any future examination 
of model spaces for Hartree-Fock theory. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, November 13, 1970, 

129-02. 
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APPENDIX - THE BAKER APPROXIMATION 


The high-energy approximation for electron scattering consists of using the ap- 
proach of Glauber to solve the Dirac equation for an electron in the presence of a scalar 
potential. Baker shows that, in this approximation, the matrix element which must be 
evaluated is 


M = V -f' J 0 (qb)[e ix ^ - l]b db 

where 

X(b) = - — / V(b, z)dz 
k *ioo 


( 1 ) 


( 2 ) 


V(b, z) = V(r) = / 3(R)V c ( |R - r |)dR (3) 

Here p(R) is the nuclear density function and V c ( |R - r |) is the screened coulomb po- 
tential, 


V (x) = - f( x ) 
x la/ 


where a is the fine structure constant, Z is the nuclear charge, and a is the screen- 
ing radius, and the screening function used here is 


f(x\ = 1 - . 

W VJ77? 


Assuming that the nuclear density is spherically symmetric, it can be shown that 


x(b) = 83/77 


r »<*> / 

•Spl -sc\ 


V Q (r, R)dz R^ dR 


(4) 


V o o 

b + z , V^r, R) is the zero order multipole of V c , and -q = »Z(E/k). 

Nov/ V Q (r, R) must be known only for values of R < Rq under the assumption that p(R) 
is negligibly small for R > Rq. Based on this assumption, it can be shown that 
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Then it follows that 


[M-D 


V 0 (r,R) =• 


r > R 


.JR" , r < R J 


x(b) = Xo( b ) + Xj(b) 


where 


Xo(b) = V 


log(l + 



and 


Xj(b) = 87 777 




i^p(R)R 2 


dR 


We have found it convenient to integrate equation (1) by parts, obtaining 


(5) 


(6a) 


(6b) 


(6c) 




e i x(b)b db 


(7) 


where x'(b) = dx(b)/db. The evaluation of xqW is simple, and it can be shown that 


x ’(b) = — 77 
b 




( 8 ) 


This can then be integrated to obtain Xj(b). Now, if we assume that p(r) can be repre- 
sented by a series of Gauss- Laguerre functions 
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( 9 ) 


p( R) = £ A n^iS a>f 3r ) 


n=0 


2_2 




and then make use of the addition theorem for Laguerre polynomials, it can be shown 
that 


i ^b 2 + r 2 ) = e" 


aj3 2 b 2 


p(r) - (/3b) 


V" L m-l^ V ) 

/. i m 

m=l 


A m+n^®> ^ r) l 


From this we obtain 


x'i(b) = ^ ?7e“ 




47T(/3b) Z > — L^ 1 _ 1 ( i 3 2 b 2 )B 

m 


ml 


where 


B 


m 


S) Am+nDn 


and 


D n ~ Jf\ ^ a ’^ r2dr 


Equation (8) may be integrated to yield 


Xl (»> = - 2 » / «'"' 32r2 ^ + 4 «) £ 

•4* R m=l 


Im(x) = X e_Q,yL m (y)dy 
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The XqO 3 ) and XjO 5 ) and their derivatives are then used in equation (7) to obtain M. 
The final evaluation of the integral is carried out using the same general procedure as 
was used by Baker. The electron scattering form factor is related to the square of M 
divided by the point cross section. The change in the electron wavelength as the elec- 
tron approaches the scattering center must be accounted for in the calculation of the 
form factor. The correction factor suggested by Ravenhall is used in our calculations 
(ref. 37). 
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TABLE L - DESCRIPTION OF BASIS SPACES TO BE USED IN 


HF CALCULATIONS 


Radial dependence 

Space 

Basis 

Harmonic oscillator 

1 

( ls l,2’ 2s l 2’ ld 3/2’ ld 5,2^ lp l 2’ lp 3'2^ 

Harmonic oscillator 

2 

( ls l, 2’ 2s l 2’ ld 3/ 2’ ld 5/2^ 

UPl 2 , lp 3 2 , 2p x 2 , 2p 3 2 , lf 5/2J lf 7 2^ 

Harmonic oscillator 

3 

( ls l 2 , 2s 12 , ld 3 , 2 . ld 5 , 2 , 3sj 2 , 2d 3 2 , 
2d 5 2’ % 2' lg 9 2^* 

( lp l 2 ' lp 3 2’ 2p l. 2* 2p 3 2’ lf 5/2’ lf 7/2^ 

Wood- Saxon 

-- 

( ls l 2 , 2s x 2 , ld 3 2 , ld 5 2 ).(lp 1 2 , lp 3 2 ) 



TABLE II. - COMPARISON OF HF WAVE FUNCTIONS. SINGLE PARTICLE ENERGIES. HF ENERGIES. 


AND GAPS FOR NEON-20 
[Oscillator length, 1. 88 fm J 


(a) HF wave functions and single particle energies 


Magnetic 

Space 




Positive parity states 




Single 

quantum 











particle 

number 


is i/ 2 

2s l 2 

3s 1/2 

ld 3 2 

2d 3 2 

ld 5 ' 2 

2d 5 2 

ig 7 2 

'* 9.2 

energy. 

m 





Negative parity states 




e A- 

and parity 











MeV 

77 

(a) 

^1/2 

2p l/2 

lp 3/2 

2p 3 2 

lf 5/2 

lf 7 ' 2 






1 

0. 9894 

0. 1336 


-0.0387 


0.0429 




-38.64 

1 

2 

2 

0. 9797 

0. 1890 


-0. 0446 


0.0490 




-44. 34 


3 

0. 9710 

0.2156 

0.0635 

-0. 0510 

-0. 0152 

0.0566 

0. 0169 

-0. 0149 

0. 0120 

-46.09 


l 

0. 3173 


-0.9483 







-19. 16 

— 









1 

2 

2 

0. 3249 

0.0110 

-0. 9265 

-0. 1454 

0.0766 

-0.0942 




-23.08 


3 

0. 3536 

0.0218 

-0.9091 

-0. 1619 

0. 0947 

-0. 1130 



— 

-24. 21 


1 




1.0000 







-15. 52 

3" 

"2 

2 



0.9752 

0.2145 

0.0332 

0.0432 




-18.43 


3 



0.9687 

0.2411 

0.0363 

0.0456 




-18.45 


l 

0. 9483 


0. 3173 







-11. 55 

— 









1 

2 

2 

0. 9315 

0. 1235 

0.3213 

0. 0879 

0.0781 

-0.0051 




-13.05 


3 

0. 9158 

0. 1429 

0. 3510 

0. 1033 

0. 0832 

-0. 0088 




-12. 72 

i + 

1 

0.0212 

-0. 5062 


-0. 3196 


0. 8008 




-4. 46 

1 

2 

2 

0. 0443 

-0. 5114 


-0. 3165 


0.7977 




-6. 06 


3 

0. 0468 

-0. 5016 

-0.0417 

-0.3577 

-0.0109 

0.7697 

0.0637 

-0.0788 

0. 1168 

-7. 96 


a For definition of spaces see table I. 


(b) HF energies and 
gaps 


Basis 

Energy, 

Gap, 


e hf- 

MeV 

A, 

MeV 

1 

-28.03 

4. 16 

2 

-36. 85 

5. 36 

3 

-40.02 

7. 03 
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TABLE in. - PROJECTED RMS RADIUS AND 


REDUCED E2 RATE 


Basis 

Oscillator 

length, 

b, 

fm 

rms radius. 

2 1 2 

(R ; 

fm 

Transition rate, 
B(E2: CT - 2 + ), 
e 2 • fm 4 


Fixed oscillator length 

1 

1.88 

2. 94 

109. 9 

2 

1.88 

2.77 

164. 5 

3 

1. 88 

2. 73 

190.0 


Fixed size of space 

3 

1. 82 

2.68 

174.7 


1. 88 

2. 73 

190.0 


1. 93 

2.78 

204. 1 


Fixed rms radius 

1 

1. 78 

a 2. 78 

82.9 

2 

1. 88 

2. 78 

164. 5 

3 

1. 93 

2. 78 

204. 1 


a Average of projected rms radii obtained with 
spaces and oscillator lengths shown. For 

2 1 2 

spaces 1 and 3, values of < R ) are 2. 79 and 
2. 77 fm, respectively. This introduces an un- 
certainty of about 2 percent in predicted E2 
rates. 


TABLE IV. - REDUCED E2 TRANSITION 
RATES WITHIN GROUND STATE 
BAND OF NEON-20 


[Space, 3; oscillator length. 
1. 93 fm] 


Transi- 

tion 

Predicted 

Experiment 

Transition rate. 
B(E2; J. - .J f ). c 2 • fm 4 

°i- 2 i 

204 

286: 15 

2 i- 4 i 

103 

89t 9 

4 i ~ 6 i 

90 

129t 13 




TABLE V. - COMPARISON OF HF WAVE FUNCTIONS. SINGLE PARTICLE ENERGIES, HF ENERGIES. 


GAPS. RMS RADII. AND E2 RATES FOR MAGNESIUM-24 
[Oscillator length, 2.09fm.J 


(a) HF wave functions and single particle energies 


Magnetic 

Space 



Positive parity states 





Single 

quantum 











particle 

number 


ls l '2 

2s l 2 

3s l 2 

ld 3 2 

2d 3 2 

ld 5 2 

2d 5 2 

ig 7 2 

lg 9 2 

energy. 

m 




Negative parity states 





e \ * 

and parity 




r- 







MeV 

77 

(a) 

lPl 2 

2 Pl 2 

lp 3 2 

2p 3 2 

lf 3 2 

lf 7 2 










-0.0554 


0.0662 




-39. 94 

1 + 









2 

3 

0. 9621 

0. 2370 

0.0672 

-0.0686 

-0.0311 

0.0819 

0.0354 

-0. 0045 

0. 0074 

-44. 14 

r 

2 

-0. 4061 

-0.0515 

0.8900 

0. 1741 

-0.0521 

0.0850 




-23.03 

2 

3 

-0. 4356 

-0. 0772 

0.8642 

0. 2037 

-0.0726 

0. 1031 




-26. 37 

3" 

2 



-0.9655 

-0. 2413 

-0.0597 

-0.0773 




-17. 41 

" 2 

3 



-0. 9550 

-0. 2754 

-0.0688 

-0.0861 




-18. 40 

r 

2 

0.8884 

0.1692 

0. 3948 

0. 1106 

0. 1064 

0.0519 




-14. 98 

2 

3 

0. 8629 

0. 2002 

0.4252 

0. 1295 

0. 1191 

0.0594 



— 

-15.69 




fl RQ97 


-0. 3024 


0.7442 




-7. 40 

i + 


U. U Do * 

-U. D 1 






2 

3 

0.0479 

-0. 5398 

-0. 1227 

-0. 3341 

-0.0323 

0.7418 

0. 1267 

-0.0511 

0.0978 

-10. 14 


2 




-0. 2823 


-0 9593 




-3. 45 

3 + 

"2 

3 




-0. 2921 

-0.0245 

v . %s U Zf 

-0.9347 

-0. 1511 

-0. 0688 

-0. 1140 

-4. 56 


a For definition of spaces see table I. 


(b) HF energies and gaps 


Basis 

Energy, 

Gap, 


E HF’ 

A, 


MeV 

MeV 

2 

-34. 49 

1. 37 

3 

-40. 87 

1. 70 


(c) rms radii and E2 rates 


Basis 

rms radius, 

2 V 2 

(R 2 ) , 

B(E2; 0+ - 2+), 
e 2 • fm 4 


fm 


2 

3.20 

351.9 

3 

3.037 

373.9 


Experiment 3.02±0.03 436*46 


TABLE VI. - REDUCED E2 TRANSITION RATES WITHIN 
GROUND STATE BAND OF MAGNESIUM -24 
[space 3; oscillator length. 2.09 fm^] 


Transi- 

Predicted 

Experiment 

tion 

Transition rate. 


P(E2; J. - 

Jj.). c ' fm 4 

°; - 2 i 

374 

436±40 

2 + — 4 + 
Z 1 

191 

133±13 





TABLE VII. - COMPARISON OF HF WAVE FUNCTIONS, SINGLE PARTICLE ENERGIES. HF ENERGIES, 


GAPS, RMS RADII, AND E2 RATES FOR SILICON-28 (OBLATE) 
[Oscillator length, 2. 09 fm.J 


(a) HF wave functions and single particle energies 


Magnetic 

Space 

Positive parity states 

Single 

quantum 

number 


ls l/2 

2 s 1/2 

3 s 1/2 

ld 3/2 

2d 3/2 

ld 5 2 

2d 5/2 

lg 7/2 

lg 9/2 

particle 

energy, 

m 

and parity 


Negative parity states 

e X’ 

MeV 












(a) 

lp l/2 

2p l/2 

lp 3/2 

2p 3/2 

lf 5/2 

lf 7/2 










i + 

2 

0. 9769 

0. 1946 


0. 0566 


-0.0678 




-42.62 






2 

3 

0. 9631 

0.2365 

0.0650 

0.0631 

0.0312 

-0.0759 

-0.0363 

-0.0109 

0.0111 

-48.66 


2 



-0.9772 

-0. 1842 

0.0595 

0.0868 




-24. 86 

3 







2 

3 



-0. 9634 

-0. 2413 

0.0679 

0.0954 




-28. 75 









2 

0. 6328 

0.0856 

0. 7461 

0. 1573 

-0.0633 

-0.0829 




-22. 99 

1“ 





2 

3 

0.6505 

0. 1335 

0.7128 

0. 1935 

-0.0755 

-0.0885 




-26. 59 







2 

0. 7482 

0. 1395 

-0.6171 

-0. 1701 

-0.0927 

0.0486 




-15. 70 

1" 





2 

3 

0.7155 

0. 1782 

-0.6350 

-0. 2039 

-0. 0955 

0. 0497 




-16. 79 






_ -t- 

2 






1.000 




-8. 78 

5 










2 

3 






0.9705 

0. 2145 

0.0606 

-0.0921 

-11. 66 








. -L 

2 

-0. 1059 

0.7928 


-0. 2375 


0. 5512 




-7. 82 

1 + 







2 

3 

-0. 1249 

0.7356 

0. 1770 

-0. 2511 

-0. 0172 

0. 5730 

0. 1069 

0.0451 

-0.0824 

-10.04 

* -I- 

2 




0. 6893 


-0. 7245 




-6. 30 










" 2 

3 




-0. 6846 

-0. 1109 

0.6950 

0. 1561 

0. 0562 

0.0922 

-8. 51 





! 


a For definition of- spaces see table I. 


(b) HF energies and gaps (c) rms radii and E2 rates 


Basis 

rms radius, 

2 1 2 
( R^ 

fm 

B(E2; Oj - 2 j). 
e 2 - fm 4 

2 

3. 336 

438. 9 

3 

3. 096 

376. 2 

Experiment 

3. 08±0. 06 

327 * 17 


Basis 

Energy. 

E HF* 

MeV 

Gap, 

A, 

MeV 

2 

-50. 30 

5. 31 

3 

-61. 16 , 

7. 24 


31 





TABLE VIII. - COMPARISON OF RMS RADIUS RADIUS AND REDUCED 
E2 RATES FOR SILICON-28 


^Oscillator length, 2.09 fm] 


Description 

rms radius, 
2 1/2 
<R > - 

fm 

2 4 

Transition rate, e • fm 


B(E2; 0| - 2j) 

B(E2; 0^ - 2^) 

Oblate 

3.096 

376. 2 

— 

Inverted co- 
existence 
calculation 

3.092 

342. 3 

34.2 

Experiment 

3. 08±0. 06 

327±17 

i 35±4 
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TABLE DC. - COMPARISON OF HF WAVE FUNCTIONS, SINGLE PARTICLE ENERGIES, HF ENERGIES, 
GAPS, PROJECTED RMS RADIUS, AND E2 RATES FOR SULFUR- 32 (PROLATE) 

[Oscillator length, 2.19 fmf] 

(a) HF wave functions and single particle energies 


Magnetic 

Space 




Positive parity states 




Single 

quantum 

number 


ls l/2 

2s 1/2 

3s 1/2 

ld 3/2 

2d 3/2 

ld 5/2 

2d 5/2 

lg 7/2 

lg 9/2 

particle 

energy, 

m 

and parity 





Negative parity states 




e \’ 

MeV 

(a) 

lPl/ 2 

2 Pl/2 

lp 3/2 

2p 3/2 

lf 5/2 

lf 7/2 




7T 







2 

0.9783 

0. 1982 


0.0372 


-0.0476 




-41. 22 

1 







2 

3 

0. 9617 

0.2540 

0.0684 

0. 0432 

0. 0240 

-0.0536 

-0.0262 

-0.0021 

0.0014 

-48,95 


2 



-0. 9773 

-0. 2038 

0.0308 

0.0485 




-24. 52 

£ 







2 

3 



-0.9581 

-0. 2786 

0.0369 

0.0554 




-29. 59 









2 

0.6130 

0. 1013 

0. 7635 

0. 1636 

-0. 0277 

-0.0590 




-22 62 

1 






2 

3 

0.6568 

0. 1706 

0.7017 

0. 2051 

-0.0386 

-0.0598 




-27. 27 







2 

0.7652 

0. 1514 

-0. 5981 

-0. 1594 

-0.0695 

0.0602 




-18.66 

1' 





2 

3 

-0. 7010 

-0. 2039 

0.6415 

0.2143 

0.0724 

-0.0658 




-21. 53 







2 






1.000 




-9.06 

5 + 










2 

3 






0.9621 

0.2665 

0.0298 

-0.0498 

-12.38 









2 

-0. 1336 

0. 8203 


-0. 1865 


0.5234 




-7.99 

1* 







2 

3 

-0. 1668 

0.7393 

0. 2238 

-0. 1911 

-0.0295 

0.5633 

0. 1327 

0.0199 

-0. 0533 

-10. 47 

„-i- 

2 




0. 5437 


-0.8392 




-7.08 

3 









- 2 

3 




0. 5557 

0. 1290 

-0.7912 

-0.2088 

-0.0146 

0.0683 

-9.64 






4 - 

2 




-0. 8392 


-0. 5437 




-4. 68 

3 + 









2 

3 




-0. 8003 

-0. 1727 

-0. 5523 

-0. 1311 

0. 0683 

0. 0535 

-6. 09 


a For definition of spaces see table I. 


(b) HF energies and gaps (c) Projected rms radius 


and E2 rates 


Basis 

Energy, 

E HF' 

MeV 

Gap, 

A, 

MeV 

. 




Basis 

rms radius, 

2 12 

'R > 

B(E2; Oj - 2p. 

2 

-51. 3 

0. 71 



fm 


3 

-64. 74 

1.02 


2 

3. 567 

401. 9 




3 (prolate) 

3. 240 

325. 1 




3 (oblate) 

3. 246 

244. 2 




Experiment 

3. 23±0. 07 

2 17±30 
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TABLE X. - COMPARISON OF HF SINGLE PARTICLE ORBITS IN NEON-20 FOR 


WOOD-SAXON (WS) AND HARMONIC OSCILLATOR (HO) BASES 


Magnetic 

quantum 

number 

Basis 

ls l 2 

2s l 2 

ld 3 . 2 

ld 5 2 

lp l 2 

lp 3- 2 

m 








and parity 









WS 

0. 9969 

0.0718 

-0.0196 

0. 0249 



1 + 



9 

HO 

0. 9881 

0. 1436 

-0. 0345 

0. 0422 







WS 





-0. 4253 

0. 9051 

1 





"o 

HO 





-0. 5774 

0. 8165 

La 





3" 

WS 






-1.0000 

"2 

HO 






-1.0000 


WS 





0. 9051 

0. 4253 

f 





2 

HO 





0. 8165 

0. 5774 







WS 

-0.0279 

0. 7009 

0. 4760 

-0. 5305 



1 + 



9 

HO 

0. 0346 

-0. 5540 

-0. 5261 

0.6443 



La 




TABLE XI. - COMPARISON OF NUCLEAR PROPERTIES WITH 
HARMONIC OSCILLATOR (HO) AND WOOD-SAXON (WS) BASIS 


[Oscillator length, 1. 924 fm.J 


Nucleus 

Basis 

rms radius. 

Transition rate. 

Energy, 

Gap. 



2 1 2 
FT 

B(E2; 0* - 2*), 

o A 

e hf- 

A. 



fm 

e ^ • fm 

MeV 

MeV 

20 No 

HO 

3.01 

135.7 

-98. 37 

8.33 


WS 

2.84 

140. 8 

-84. 98 

8.04 

2 V 

HO 

3. 12 

191. 2 

-116.77 

0.24 


WS 

2. 96 

174. 4 

-87. 54 

0. 11 

ro 

CO 

W 

HO 

3. 20 

269.8 

-160. 15 

9. 11 


WS 

3.09 

274. 9 

-116. 03 

8.62 

32 s 

HO 

3. 25 

192. 8 

-193. 75 

0.24 


WS 

3. 17 

182.7 

-134. 15 

0. 22 


Single 

particle 

energy, 

e A’ 

MeV 

-57. 95 
-47.69 

-28. 10 
-25. 13 

-23. 39 
-19. 17 

-21. 63 
-19. 17 

- 6 . 10 
-9. 97 
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Oscillator length, b. frc 

igure 1. - Plot of HF energy as function of oscillator length for 
neon-20 (space, 2). 



Figure 2. - K eon-20 elastic form factor for 200-MeV electrons. Oscillator 
length, 1. 93 femtometers; basis space, 3. 




Figure 3. - Neon-20 elastic form factor for 200-MeV electrons. Comparison 
for 6, 10, and 15 state spaces. Proton radius, 0. 77femtometers. 


Scalar form factor, f, 



Figure 4. - Neon -20 scalar form factor. (LSJ = 202. ) 




O Experiment (unpublished data 
from R. de Swiniarski) 

Glendenning-Veneroni 

Thompson-Tang 

O O 

O 



40 60 80 100 120 140 160 


Center of mass angle, deg 

Frure 7. - Predicted angular distribution for inelastic proton scattering com- 
pared wit> experiment "using form factors from neon-20 from figure 6. Pro- 
ton lab encv f v , 24. 5 iV.eV', Q -value, -1. 63 IV.eV (2 + l 


Energy, E HF 


-36 



1.72 1.76 1.80 1.84 1.88 1.92 1.96 2. 

Oscillator length, b, fm 


Figure 8. - Plot of HF energy as function of oscillator length 
for magnesium-24. Ten state space. 





Center of mass angle, Oq^, deg 

Figure 10. - Predicted angular distributions for inelastic proton scattering 
from magnesium -24 compared with experiment for 15 state space. Oscil- 
lator length, 2. 09 femtometers; proton lab energy, 17. 5MeV; Q-values, 
1. 37 MeV (2 + ). 




O Experiment iref. 34) 



$ 


<5 



Center of mass angle, 0 CM , deg 

Figure 12. - Predicted angular distribution for inelastic proton scattering from 
silicon-28 compared with experiment Oscillator length, 2. 09 femtometers; 
basis space, 3; proton lab energy, 17. 5MeV; Q -value, -1. 78MeV (2 + ). 


Form factor, |F(q)| 



Figure 13. - Sulfur-32 elastic form factor for 262 -MeV electrons. Oscillator 
length, 2. 49 femtometers; basis space, 3. 



O Experiment (ref. 34) 


-5 

5 « 

4 ? 



•° 30 50 70 90 no 130 13 

Center of mass angle, 0 CM , deg 

Figure 14 - Predicted angular distribution for inelastic proton scattering from 
sulfur-32 compared with experiment. Oscillator length, 2. 19 femtometers- 
basis space, 3 : proton lab energy, 17. 5 MeV; Q-value, -2. 24 MeV (2 + l 




fil 



Figure 16. - Comparison of elastic form factors for sulfur-32 using Wood-Saxon 
and harmonic oscillator bases with 6 states. Energy, 262 MeV. 


Cross section, mb/sr 


20 r- 


O Experiment (unpublished 

data from R. de Swiniarskil 



20 40 6 0 80 100 120 140 160 

Center of mass angle, 0 CM , deg 


Figure 17. - Predicted angular distributions for inelastic proton scattering from 
neon-20 compared with experiment using Wood-Saxon and harmonic oscilla- 
tor bases with 6 states. Proton lab energy, 24. 5 MeV; Q -value, -1. 63 MeV (2 + ). 
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